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Abstract. The photon plus heavy-flavour quark (6, c) final state provides a unique and 
valuable window into both the sea quark content of the proton and the splitting of gluons 
into heavy-flavour quark pairs. A new combination of experimental techniques has provided 
the basis for the first measurements of the differential 7 -I- c+ X and 7 -I- 6-1- X production cross 
sections at = 1.96 TeV. The measurements use ~1 fb~^ of data from pp collisions collected 
with the D0 detector. 



1. Introduction 

The results presented are the first measurements of the differential 7 -|- c jet and 7 -|- 6 jet 
production cross sections at a hadron-hadron collider [1]. The measurements used ~ 1 fb~^ 
of data and were performed at the Fermilab Tevatron pp collider at ^/s = 1.96 TeV collected 
using the D0 detector. The presented production cross sections are differential with respect 
to the photon candidate's transverse momentum (pj^) and are divided into five pj. bins. The 
measurements are also differential with respect to the leading jet's rapidity (y-''^*) and the leading 
photon's rapidity (y^) and are binned in two regions of photon- jet rapidities: y'^ ■ yJ^* > and 
yiet . y7 < 0. 

2. Data Selection 

Events are recorded using the D0 detector [2] and must contain at least one photon candidate 
and at least one heavy-flavour jet candidate. The photon and jet candidates with the highest 
transverse momenta (leading) are selected. The pT of the photon candidate must be greater 
than 30 GeV/c with \y'^\ < 1.0. The leading jet must satisfy pr > 15 GeV/c and \y^'^^\ < 0.8. 
The primary collision vertex must be located within 35 cm of the detector's center along the 
beam pipe, reconstructed with at least three associated tracks. To suppress background events 
coming from cosmic-ray muons and W ^ Iv decays, the total missing transverse energy [E^^'^'^''^) 
in the event must satisfy the condition of i?-^^**** < 0.7 • pj.. 

This analysis employs triggers that distinguish large depositions of energy in the 
electromagnetic (EM) calorimeter, and events are selected if at least one of these triggers has 
fired. The photon candidate is reconstructed, forming a cluster of energy. The reconstruction 
procedure sums the total energy deposited in the calorimeter within a cone of radius TZ = 0.2, 
where TZ = ^ (Ar/)^ -|- (A(/))2. At least 96% of the photon candidate's energy must be in the 
EM calorimeter, and it must be spatially well isolated. To be considered isolated, the total 
energy in a cone of 7^ = 0.4, excluding the candidate's energy in the EM calorimeter (-Eem), 
must be less than 7% of -Eem- Additionally, the energy deposition within the third layer of the 



EM calorimeter must be consistent with that of an electromagnetic shower. Photons traversing 
the D0 detector typically do not leave tracks, thus the photon candidate is required to have 
less than a 0.1% probability of being spatially matched to any track in the event. To reject 
additional backgrounds coming from dijet events, an an artificial neural network (7— ANN) [3] 
is used with the requirement of 7— ANN > 0.7. 

Jets are reconstructed using D0's Runll algorithm [4] with a radius 0.5. The leading jet is 
subjected to the additional requirement of having at least two associated tracks with hits in the 
silicon tracker, one with pT > 1.0 GeV/c and the another with > 0.5 GeV/c. Heavy-flavour 
jets are identified using information arising from the long lifetimes of their hadrons. To best 
enhance the fraction of jets originating from heavy-flavour quarks, an artificial neural network 
(6— ANN) is employed. The 6— ANN [5j is trained such that the output of heavy-flavour jets 
tends toward one and that of light {u, d, s quarks and gluons) jets tends to zero. For this 
selection, the requirement of 6— ANN > 0.85 must be satisfied. 

3. Purity Estimates 

The fraction of direct photons in the final data sample is determined from a simultaneous fit of 
the 7— ANN output shape in data using templates derived from Monte Carlo (MC) signal photon 
and background dijet events simulated using pythia [6]. The fitting procedure determines the 
linear combination of the signal and background templates that best describes the output shape 
in data. The fit is performed in each pj, bin for both rapidity regions separately. 

The fraction of c and b jets in the final data sample is determined using a similar method 
to that of photon purity technique. However, for the heavy-flavour fraction fitting technique, 
another discriminant, Pb-jct was used. Pb-jet is defined as -Pfe-jet = ~ liillj ^tmck' where -Ptrack 
the probability of a track in the jet cone to originate from the primary collision vertex, omitting 
the least likely track to have come from this vertex. The value of Pb-jet is large for b jets, and 
tends toward zero for light jets. The shape of the -Pfe-jet distribution in data is fit from MC 
templates for both b jets and c jets, with the light jet template coming from a light jet enriched 
sample in data. The fit is performed in each p^ bin separately, and the quality of the fit in each 
bin is of good quality. As an example, the result of the fit for 50 < < 70 GeV/c is shown in 

Fig.m 

4. Cross Section Results 

The measured differential cross sections are presented in five pj, bins and two regions of y'^ ■ y-''^*, 
and can be seen in Fig. [2] for both the 7 -|- 6 and ■y + c processes. Due to the finite resolution of 
the calorimeter, pj. smearing effects were accounted for using the unfolding method described 
in Ref. [7]. Statistical uncertainties for these results vary between 0.2% for 30 < pj, < 40 GeV/c 
to ~9% for 90 < pj, < 150 GeV/c, and the systematic uncertainties range from 15% to 28%. 
The theoretical predictions from next-to-leading (NLO) order calculations |8] are also presented 
in Fig. [21 where the renormalization scale ^/j, factorization scale fip and fragmentation scale fif 
are all set equal to p'^. 

The ratios of the measured to the predicted cross sections for both j + b and 7 -|- c cross 
sections in the two rapidity regions are shown in Fig. [3l The uncertainty due to the scale 
choice is estimated by simultaneously varying all three scales by a factor of two, Ai_R,F,/ = 0-5p^ 
and IJ'R^FJ = "^Pt- The CTEQ6.6M parton distribution functions (PDFs) were used in the 
theoretical calculations, and their corresponding uncertainties were calculated according to the 
prescription in Ref. [9]. Agreement within uncertainties can be seen in Fig. [3] for the 7 -|- 6 
cross sections between data and theory in both rapidity regions. The 7 + c cross sections show 
agreement within uncertainties up to pip ~ 70 GeV/c; however, there is a rising discrepancy 
between data and theory that grows as pip increases. Certain models, including those that 
include an intrinsic charm (IC) component to the proton, predict higher cross sections for 7 -|- c 
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Figure 1. Distribution of the Pb-jet discriminant in the final data sample in the range of 
50 < pj^ < 70 GeV/c. The templates for light, c and b jets have been normalized to the number 
of events in data and then weighted by their found flavour fractions, and the shape of their 
sum is found to be in good agreement with that of the data. The shown uncertainties for the 
data are statistical in nature, whereas for the templates they incorporate statistical and fitting 
uncertainties. 



> 

Q. 



■o 

Q. 



10 
1 

10-1 



D0, L. = I.Ofb"^ 

int 



Y+c+X 



S 10' 



10" 



10' 



^ Y+b+X 



|y|<i.o 
l/'l < 0.8 
p*^^' > 15 GeV 





_L 



• yV" > 
■ yV" < 

NLO QCD 
CTEQ 6.6M 

U, = p' 

R,F,f T 



(x3.0) 
(x1.0) 

(xO.3) 

(xO.1) 





20 40 60 80 100 120 140 



p; (GeV) 



Figure 2. The 7 + 6 and 7 + c differential cross sections as a function of pj^ for both rapidity 
regions. The data points include the overall uncertainties from the measurement, and the 
theoretical predictions are displayed as dotted lines. 
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Figure 3. The ratios of the measured to the predicted cross sections for 7 + 6 production (top) 
and 7 + c production (bottom) in each rapidity region. The data points include both their 
statistical (inner tick) and overall uncertainties (entire error bar). The uncertainties from the 
theoretical predictions include those from the CTEQ 6.6M PDFs (yellow band) and from the 
choice of scale (full line). The ratio of two intrinsic charm models to the standard theoretical 
predictions are also included (dashed lines). 



production. Although at least one IC model shows better agreement with the data than the 
standard theoretical predictions do, neither model shows agreement for the entire pj. range in 
either rapidity region. 
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